Quantum polarization spectroscopy of ultracold spinor gases 
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We propose a method for the detection of ground state quantum phases of spinor gases through 
a series of two quantum nondemolition measurements performed by sending off-resonant, polarized 
light pulses through the gas. Signatures of various mean-field as well as strongly-correlated phases 
of F = 1 and F = 2 spinor gases obtained by detecting quantum fluctuations and mean values of 
polarization of transmitted light are identified. 
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It has been demonstrated several years ago that funda- 
mental quantum spin noise of a collection of cold atoms 
can be measured via quantum noise limited polarization 
spectroscopy [l[ . Since then, a quantum interface of light 
with optically thick atomic spin ensembles has become a 
promising and powerful method for a transfer of quan- 
tum information between atomic internal degrees of free- 
dom and light. The basic concept underlying such atom- 
light interfaces is provided by the off-resonant coupling of 
the collective atomic spin, i.e., of several magnetic sub- 
levels, to the polarization of light. In particular, such an 
off-resonant interaction, followed by a quantum measure- 
ment on light, has been shown to be a powerful quantum 
nondemolition (QND) tool to generate spin squeezed and 
entangled atomic states ill IS , to tclcport quan- 
tum states between ensembles [J, or to propose [H 0] 
and realize Q high fidelity quantum memories for light. 



In the present paper we propose to apply quantum 
polarization spectroscopy techniques for the detection of 
various quantum phases of degenerate atomic gases, i.e., 
with atoms having spin degrees of freedom. Such ultra- 
cold spinor gases have recently brought a new perspective 
to the study of magnetic systems. Seminal experiments 
of the MIT group with an optically trapped spin F = 1 
[lot , a nd theory papers of Ho 



Sodium condensate [lOl ] , a nd theory papers of Ho [llj . 
and Ohmi and Machida [l2j have triggered the use of cold 
atoms to study magnetic ordering and domains. Spin in- 
teraction effects are very much enhanced in the strongly 
correlated regime, which nowadays is reachable experi- 
mentally [l3l | by loading an ultracold spinor gas into an 
optical lattice so that the kinetic energy (tunneling) be- 
comes small in comparison with atom-atom interactions. 
Then the atoms can be well described by a generalized 
spinor Bose-Hubbard Hamiltonian (BHH) 0, [U, 
In the limit of small occupation number it reproduces 
accurately (within the experimentally achievable regime) 
some of the most paradigmatic spin chain models. Ex- 
perimental observation of the rich variety of magnetic 



ordering present in these systems remains however elu- 
sive due to similiar values of the scattering lengths on the 
different spin collision channels. 

A way to determine properties of a quantum spinor gas 
could be to use a strong QND measurement. As shown 
here, a series of QND measurements using polarization 
spectroscopy on light transmitted through a condensate 
yields mean values and variances of the atomic total spin 
operators, thus allowing unambiguous distinction of var- 
ious atomic quantum phases. 

Formalism - We consider a sample of neutral atoms in 
a 2F+1- dimensional ground state manifold \F, m) , inter- 
acting off-resonantly with linear ly p olarized light propa- 
gating along the ^-direction (cf. [17j). After adiabatically 
eliminating excited atomic states, the interaction can be 
described via an effective Hamiltonian 
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where L is the length of the atomic sample. The condi- 
tions under which decohcrence due to absorption of light 
can be neglected, such that this Hamiltonian is valid, will 
be discussed below. In Eq. (Q]), ai oc h'y\ 2 c/(16irAA) 7 
where A is the cross-section of the atomic sample over- 
lapping with the probe light, A is the detuning, p is the 
(in general ^-dependent) atomic density, A is the wave- 
length, and 7 is the excited state line width. s a = s a (z, t) 
(s± = s x ± isy) are the components of the Stokes vector 
characterizing the polarization of the light pulse, 4>(z,t) 
is the photonic density, and j = j(z,t) are atomic spin 
operators. The term proportional to ao corresponds to 
the AC Stark shift, while a 2 — > for values of A large 
compared to the excited state hyperfine structure [171 ]. 
Here we assume a 2 = and restrict to the linear cou- 
pling between the Stokes operator and the atomic spin, 
which represents a QND Hamiltonian. Using Heisen- 
berg equations of motion, we find that j z is conserved. 
For a pulsed probe with a duration of fis we can thus 
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assume j z (z,t) = j z (z), as spin diffusion happens on 
a much larger timescale of ms. For a probe strongly 
polarized along the x-direction, the macroscopic Stokes 
operator S x (S a = J dts a ) can be replaced by a c- 
number S x sa (S x ) = Np/2, being Np the number of 
photons. Neglecting retardation effects, the propaga- 
tion equation for S y reads d z S y (z) = —aiAp(S x )j z (z). 
It is convenient to introduce the collective spin in z- 
direction J z = L dz Apj z [z) and to define quadrature 
operators via X$ = \f2/NpS yi Yg = y/2/NpS z , such 
that [JfgjYg] sa i. Integrating the propagation equation 
gives 



X" out 
S 



^FWa 



(2) 



with k 2 = a 2 NpNA^ (number of atoms Na = J dz pA). 
Eq. © is valid provided that 2n 2 FN A /N P < 1. Fluctu- 
ations of the Xs quadrature are calculated as 
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where 1/2 stems from the quantum fluctuations of the 
coherent input state of light (preparing a squeezed input 
state would allow to reduce this contribution). 

Equations |TH3]) allow to estimate the feasibility of the 
proposed methods for measurement of the fluctuations of 
J z of a spinor condensate. A strong QND measurement 
repeated on samples prepared in the same state yields a 
complete knowledge about the operator J z . The strength 
of the interaction is determined by the constant n. In or- 
der to obtain useful information about the spin fluctua- 
tions, the second term in Eq. ([3]) must be large compared 
to 1/2, the quantum noise of a coherent probe. For the 
sake of this estimate we take a ferromagnetic spinor con- 
densate (also referred to as a coherent spin state outside 
the context of quantum gases), for which the second term 
is equal to k 2 /2 [H, 0]. The interaction constant can be 
expressed [l8j as k 2 = arj, where rj is the probability of 
spontaneous excitation caused by the off-resonant probe 
and a is the resonant optical depth of the sample, i] 
describes the decoherence which must be kept small by 
choosing the detuning and/or the strength of the probe 
pulse, in order to minimize the distortion of J z and to val- 
idate the use of Hamiltonian Eq. (fT]). The optimal value 
for a QND measurement of J z is n 2 = y/a/2 18 1. Hence 
for an optically thick spinor condensate with a 1 the 
proposed method works well. For degenerate gases, typ- 
ical values of a = 300 or even more have been demon- 
strated [lj|. 

The most general Hamiltonian describing a spinor 



atomic sample is given by 
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+ U trap (r) * m (r) + Vint, 



(4) 

where the field operator ^^(r) creates a particle with 
spin projection m at position r and Vint describes the 
atom-atom interactions. For two identical spin-F bosons 
interacting via s-wave collisions, VS n t = Em=o9 m ^ 2B ' 
with P2M being the projector onto the subspace with to- 
tal spin F tot = 2M, and g2M the interaction strength 
for the given spin channel. For a spatially uniform trap- 
ping potential, condensation occurs in the zero momen- 
tum state, and variation al g round states are obtained 
from a trial wavefunction 12011 



10 = 



1 



F 

E ' 

n=-F 



N A 



vac) 



(5) 



with the complex components of the vector £ as param- 
eters, at creates a particle with spin projection m 
in the k = state. Ground states have been discussed in 
detail for F = 1 [ll|, F = 2 (2^1, and F = 3 Given 



the spinor £, means and variances are calculated as 
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In the strongly correlated regime, Eq. ((4| reduces to a 
spinor BHH. For an integer filling factor and considering 
tunneling pcrturbatively up to second order, an effective 
Hamiltonian with nearest neighbor spin-spin interaction 
arises. For such a situation, ground state configurations 
have been analyzed in [3] for F = 1 and in la, 22 1 for 
F = 2. Eq. (J3j) reduces then to a sum over correlations 
between all pairs of atoms: 
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Notice that in H^ t , Eq. (fT]), light couples to individual 
atoms. Thus for more than one atom per site, j k denotes 
the spin operators of the kth atom, not the total spin at 
site k. 

The procedure described so far allows to obtain mean 
and variance of J z . Reading an orthogonal component 
of the atomic spin vector allows for a better discrimina- 
tion of distinct phases. Since the interaction with the 
first light pulse modifies the atomic spin state, here we 
assume to have at our disposal a second, identically pre- 
pared atomic sample. A lightpulse incident in the y- 
direction, again strongly x-polarized, corresponds to ex- 
change j z — > j y , jy — > — j z in Eqs. (TTH3]). Thus it allows 
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Figure 1: (color online) Possible combinations of additional fluctuations e| and efj imprinted on the light for the ground state 
phases of the F = 1 spinor gas in a uniform trap and in an optical lattice (a) and for F = 2 atoms in a uniform trap (b) and in 
an optical lattice (c). Filled areas denote cases where the mean of {X$ ut ) and/or (X[} ut ) is (generically) non-zero. The spheres 
in (a) illustrate the directions of the spinor for the extremal points of the ferromagnetic phase. 

(ia) Ferromagnetic m = 2 - £ = {e^s^, 2e~ 1 - 0! s%C[ } , ^\Sl,2e ici s c%,e ia c%) 
(Xg ut ) = -V2NCp, <X& ut ) = -V2NS a Sp, e 2 s = i« 2 Sg, e 2 v = \£ (l - S 2 a S 2 g ) 

(ib) Ferromagnetic m = 1 £ = (2 e - ia s|^, -e- ia s 2 g {l + 2C »> ~\fli S f> C l3> -e _a 4(-l + 2Cg), 2e 2ia sgc 3 g ) 

I youtx _ VW n I x>°ut\ _ vTV q C ^2 _ 5 „2 q2 2 _ 5 „2 /i q2 q2\ 

(ii) Polar £ = (e^SgS^, e^S e C in V2Cg, -e^SgC^, e^SgS^) 

{Xs ut } — 0, {Xu ut } = 0, e| — ^K 2 Sg(l + 3S 2 ), e 2 / = -|k 2 (3 — 2.S 2 + 3C 2 C|,C| — V3C r pS2gS^) 

(iii) Cyclic - as [20], Eq. (60-62), for zero magnetic field: (Xg ut ) = 0, (Xg- Ut ) = 0, e| = k 2 , e 2 / = k 2 

Table I: Spinors £ for the mean field ground state phases in the spin-2 case [2(j and corresponding results for means {Xg/ij} 
and additional fluctuations es/u- Abbreviations s x = sin |, c x — cos |, 5^ = sin a;, C x = cos a; are used. 



to read out the y-component of the atomic spin vector, 
with the corresponding quadrature denoted as Xjj. 

Detecting spin-1 quantum phases - In the mean 
field case, a gas of F = 1 atoms has two possi- 
ble ground states [ll[: a ferromagnetic one with all 
spins having maximal spin projection in some direc- 
tion (cos a sin /3, sin a sin /3, cos /3), and thus (j) 7^ 0, and 
a polar one with (j) = 0. Inserting the ferromag- 
netic spinor £ f = (e~ JQ cos 2 f, v^cos | sin |, e 4Q sin 2 |) 
into Eqs. © leads to (X% ut ) = -kVAUcos/3, 
(I™') = K^TVlsinasin/?, (AX° ut ) 2 = (1 + k 2 sin 2 /3)/2, 
(AX^) 2 = (1 + k 2 [1 - sin 2 /3sin 2 a])/2. These equa- 
tions arc valid provided the mean polarization of the 
probe light remains to be x-polarized which is true 
under the feasible assumption Na <C Np. For 
the polar phase, the spinor can be parameterized as 
£ p = (— e~~ ia sin/3, \[2 cos/3, e la sin/?)/-\/2, and we have 
(X| ut ) = = (A^} ut ) and (AA§ ut ) 2 = (1 + 2k 2 sin 2 /3)/2, 
(AX£ ut ) 2 = (l + 2K 2 [l-sin 2 asin 2 /3])/2. Characterizing 
the atomic phases by the additional noise (AX^jj) 2 — 
i = <? s iu imprinted on the X quadratures, we obtain 
for the ferromagnetic phase e 2 / = |k 2 — e|sin 2 o;, with 



< e| < ^K 2 , and for the polar phase efj = k 2 — e| sin 2 a, 
with < e| < n 2 . Possible values of the additional 
noise lie in non-overlapping triangles in the (e|, em- 
plane, see Fig. Q] (a). Thus both phases can be distin- 
guished through the noise imprinted on the light. In this 
particular case, this is also possible by comparing the 
mean values (Xu/s) (f° r i a ^P) 7^ (0, 7r/2)). 

For the F = 1 lattice gas with a single particle per 
site, the effective Hamiltonian is H\ at = Yliki) (cos <5 j fc j' + 
sin 5 (j k J 1 ) 2 ), where the sum runs over nearest neighbors. 
We consider only t he g round states of 23 Na discussed by 
Imambekov et al. [14| (for numerical investigations see 
[2^): (i) for a fully polarized state the properties of the 
out-going light are as in the ferromagnetic case discussed 
before; (ii) a state mixing total spin F tot = and F tot = 
2 on each bond, constructed as Yi m l£p)m, gives results 
as for the polar mean field state; (iii) singlets (dimers) 
can be put on every second bond (in ID), breaking trans- 
lational symmetry. The reduced on-site density matrix is 
p k = 1/3, and thus (Xg ut ) = 0. As (f k jf) = § for k = I, 
— I for nearest-neighbors in a singlet state, and other- 
wise, for an even number of sites that the fluctuations are 
unchanged: ( AXg ut ) 2 = 1/2. For an odd number of sites 
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in ID, or randomly oriented dimers in 2D, n atoms will 
be unpaired and thus (AXg ut ) 2 = 1/2 + 2n K 2 /(3N A ). 
Due to the rotational symmetry of the singlet, the same 
result is obtained for the Xjj quadrature, (iv) In ID and 
for S = — arctan(l/3), the ground state is a valence bond 
solid (VBS) state [24|. In this case, two-site correlations 



decay as (j k J z ) = (* + I) UM- Since the VBS 



is non-magnetized, the means of the quadrature compo- 
nents remain unchanged. As the sum of (j z j l z ) over all 
pairs of atoms gives 2/3 independent of the number of 
sites, there is no detectable change in the fluctuations: 
(AXg ut ) 2 = | + 2k 2 /3N a . Thus distinguishing between 
a dimer and a VBS state is difficult with this method, 
but in principle possible for small N A . 

For two atoms per lattice site, in the limit of vanish- 
ing tunneling the ground state consists of non-interacting 
singlets on each site. As tunneling is increased, on-site 
states with total spin 2 become important. Using a varia- 
tional ansatz ]J k \tp) k with \ip) k = cos9|F tot = 0, F* ot = 
0)+sine|_F tot = 2,F* ot = 0), asharpjump of sin from 
to a non-zero value is found as tunneling is increased 
. Evaluating the quadrature operators of the outgo- 
ing light, we find (i° s ut ) = = (X° ut ), but modified 



noise properties: (AVg ut ) 2 = (1 + 3k 2 sin 2 6 sin 2 j3)/2, 
(AX£ ut ) 2 = (1 + 3K 2 sin 2 e[l-sin 2 asin 2 /3])/2, where 
a, (3 parametrize the direction of the component with 
F tot = 2. The sharp change in the nature of the ground 
state manifests clearly in the noise properties of the out- 
going light. Let us emphasize that, as light couples to sin- 
gle atoms, here the fluctuations are different from those 
arising from fully polarized F = 2 atoms, as will be dis- 
cussed now. 

Detecting spin-2 quantum phases - For F — 2, there 
are three different ground state phases in the mean field 
case [2^: (i) ferromagnetic, (j) ^ 0, characterized by 
spin projection m = 2 (a) or m = 1 (b) in some direction, 
(ii) polar, characterized by a = S m (-l) m £m£-m 7^ 
and, (j) = 0, and (iii) cyclic, having a = and (j) = 0. 
Our results are summarized in table U and Fig. [ljb). 
Ferromagnetic and polar phases can be perfectly distin- 
guished while the cyclic phase lies within the polar one. 

Finally, we consider the F = 2 lattice gas with 
a single particle per site, following the discussion in 
[lit . The effective Hamiltonian contains up to the 
fourth power of the Heisenberg interaction: H\ at = 
J2{ki) Ea=i ^a{j k 3 l ) a - We will discuss here only some of 
the possible ground state phases arising for various com- 
binations of A Q : (i) for products of cither ferromagnetic 
or polar on-site states results are as in the mean field case; 
(ii) an anti-ferromagnetic state |£ m , £_ m , £ m , . . .), with 
|£ TO ) (|£- TO )) being the on-site state with spin projection 
m (—m) in some direction. As such Neil states are non- 
magnetized, the means of the outgoing X quadratures 
remain unchanged, but fluctuations are as for ferromag- 
netic states. Neel states can thus be distinguished from 



ferromagnetic ones due to different means, but cannot 
always be distinguished from polar states, (iii) certain 
combinations of X a favor product states |£, £, . . .) fulfill- 
ing (singlet |£, £) = 0, corresponding to a = 0. The cyclic 
states discussed in the mean field case are particular in- 
stances of this phase, which however is larger because 
the spin projection is not restricted to zero. Possible 
combinations of fluctuations for the different phases are 
summarized in Fig. [1] (c). 

Conclusions - We have presented a method to de- 
tect different magnetic orders arising in ultracold spinor 
gases, it is based on an atom-light interface, realized 
through the off-resonant interaction of a strongly polar- 
ized pulse of light and an ultracold ensemble. Fluctua- 
tions imprinted on the quadratures of the outgoing light 
contain information on the total atomic spin components. 
By comparing the fluctuations arising when the light im- 
pinges on the atomic sample from different directions, 
most ground state phases can be distinguished. We plan 
to apply our analysis to the complete classification of the 
mean field states for arbitrary F presented recently by 
Barnett et al. (22[. 
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